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Abstract 

A simple reaction-rate model is proposed for turbulent diffusion flames. 
This model was derived from the results of Chung (1972), and has reaction rate 
proportional to turbulence mixing rate. The reaction rate is also dependent on 
the mean mass fraction and the mean-square fluctmtion cf mass fraction of 
each reactant. Calculations are compared with experimental data of Kent and 
Bilger(l973) > ancL are generally successful in predicting the measured 
quantities . 
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Homenclature 

Stoichiometric coefficients of fuel, oxygen, and product, 
respectively 

mass fraction of specie ij c^, turbulent fluctuation 

empirical constat 3, equal to 1.43 and 1.92, respectively 

enpirical constants, equal to 3*0 and 0.20, respectively 

empirical constant, equal to 0,09 

diameter of Ijydrogen injector 

mean-square fluctuation of mass fraction, c^ 

static enthalpy 
stagnation enthalpy 

mean turbulence kinetic energy, u^u^/2 
press'ure 

radial distance from the axis of sjTometry 
velocity in x, y dire . bions 
turbulence velocity in x, y directions 
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U det mass-'average velocity 

J 

mole cellar weight of specie i 

reaction rate of specie i 

X, y Cartesian coordinates 

e dissipation rate of tiiThtilence 

eddy viscosity 

p density 

turbulent Prandtl-Schmldt number for transport of quantity j; 
"f “ “g " 'h ' ''t = 1-0; "c = i-3 

(“) mean or time -average value of ( ) 


Subscripts 

c 

e 

f 

J 

N 

o 


centerline value 

free -stream conditions outside mixing region 

V 

fuel 

fuel jet at injection location 

nitrogen 

oxygen 


Introduction 

There have been a nunber of studies directed toward modeling the 
effect of turbulence on chemical reactions. These studies have been motivated 
by the need to develop methods for computing turbulent, chemically-reacting 
flew fields for a wide variety of applications. Spalding proposed three such 

1--5 

models , all for cases where the chemical reaction time is much less than 
the characteristic time for turbulent mixing. His first model^ was for 
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initially untaixed reactants that react instantly on contact. It was assumed 
that the instantaneous fuel-air ratio at a point was equal to the mean value, 
plus or mi:nis the root -mean-square fluctuation. For each of these two values, 
the temperature and canposition for complete reaction was con^ted. 

Calculations indicated a finite flame width, which is qualitatively correct. 

This model has not been extensively conpared with esqperhoent, however. Spalding's 
second and third models ’ were developed for pr^ixed flames, and were 
called eddy breakup models. Each had the mean fuel reaction rate propozrtional 
to e/k, the characteristic mixing rate of the tuzbulence. The second and 
third models also had mean fuel reaction rate proportional to mean fuel mass 
fraction and to root -mean-square fluctmtion of fuel mass fraction, respectively. 
Both of these models have had some s\iccess vhen compared with escperimental 
data for premixed flows. 

It 

-^Shodes, Harsba and Peters did not use species - conservation eqxiatioias, 
and so did not model the reaction-rate terms directly. They divided the flow 
field into classes and zones, characterized by Instantaneous and mean concen- 
trations of the fuel, respectively. An assumed probability density for concen- 
tration fluctuations was used to develop a model for their effect on mean 
density. (Correlations involving density fluctuations were neglected, however). 
This model was conbined with both equilibrium chemistry and finite-rate 
chemistzy models. In the latter, each class was assumed to behave as a one- 
dimensional, transient stirred reactor. Their cooptations were compared With 
the data of Kent and Bilger^ for U-/u “10. 

J G 

6 

Hilst, et al , developed second-order cloBure methods for the species 

continuity equations, and combined these with; turbulent diffusion models. 

These models conformed to the principles of invariant modeling. A limited 

number of calculations wei*e perfom«d for two-dimensional diffusion and . 

7 

isothermal reaction of two species . Libby considered the equations 
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describing turbulent reacting flows with fast chemical react .'ons. Effects of 
density fluctuations were included in the formulation. Specific closure models 
were not .proposed, however. 

There have also been a nutober of studies which considered sii^plifled 

flow fields, in an effort to describe the interaction between turbulence and 

8 9 

nhpm-ioAi xeactions in greater detail. Among these are Chung , Alber and Batt , 
and O'Brien^*^. Ih these analyses, the probability density functions for 
velocity and scalar fluctuations were modeled and used to compute various 

U 

8 

'correlations. The analysis of Chung^ formed the basis of the model described 
here. 

This paper two purposes. First to propose a sisple reaction^rate 
model applicable to turbulent diffusion flames,’ for which the time scale for 
reaction is much smaller than the time scale for toizbulent mixing. 
Second, to provide a* limited assessment of the proposed model by comparing 
calcuIntiotU with estperimental data. 

Reaction-Rate Model 

g 

In this section, Chung's analysis of homologous shear flow forms the 
basis for a model of the chemical reaction rate in a turbulent flow. It is 
ass-limed that the chemical reaction time is much less than the turbulent mixing 

time. Also, it is assumed that the Reynolds number is large and that 

\\ 

\\ ' J 

solid boundaries have liegligible influence, so that molecu lar transport can be 
neglected. 

O 

Chung° considered a one-dimensional shear flow with nonzero derivatives 
in the y-direction, only. He also ass'imed constant density, so that the 
velocity field was not affected by the tenmerature field. For these conditions. 
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the mean species conservation equation for the ^tlel *(f) is: 

(T!^) = ^ (1) 

8 

In his solution, Chung represented the probability density o" the fluid 
elements by the sum of two half ^iaaewellian functions containing a total of 
four parameters. A nioment method was used 'solve for these parameters. 

These in turn determined mean velocity gradient, turbulence energy, and 
Reynolds stress . 

Kext, consider that the combustion is described by the following one- 
step chemical reaction: 

a (fuel) + b (oriidizer) “♦ d (product) (2) 

'I. 

i', 

■ -i 8 

Chung's solution for the probability densities, when the chemical reaction 
time is much smaller than the tvirhttlent mixing time, allows the correlation 
vc^ to he evaluated. It is: - . 

r *1^*5 if (i 

VC^ « [(1 - TT/8)k/3TTj Cj ' (3)JI 

I ■ I 

\\ 

In the combustion region, the solution for C^ is linear in y. Although , vi 
ChTong did not escplicitly d|o so, an equation for w^ can he obtained by 
substituting (3) and the solution for C^ ir.to (:^ ^ After some manipiCation, 
the resxilt can be e 3 q>ressed as : 

^ - (3/8n) 2°*5 [c^ + + (W^a/M^b) (C^ + W 

where -t is the local integral scale of the turbulence. Ih (4), the rate 
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is coMistenti with the coxjdition that tiirhulent motion controls the 
reaction rate. !Hie quantity in brackets provides an estimate of the effect of 
the pjmount of each reactant on the reaction rate* 

It is here hypothesized that (4) is approximately valid at each point 
in more ccxaplex flows than the one -dimensional shear flow considered by 

g 

Chung . This hypothesis was incorporated into a calculation method that uses 

12 

differential equations for k and e to describe turbulence transport , 
if the estimate e = is made (4) can be written in terms of 

k and < as: 

f; 

^ - A(«A)‘ [Of + * (WjA/W^b) > So°'^)] (5) 

where A is about 0.17* Comparisons with experimental data^ to be discussed 
in the next section^ suggest that A - 0.22 is a better estimatr;. 

The dependence of w^ on mass fraction resenibles a result derived by Gibson 
and Libby for steady flow relative to a reaction zone, with oxidizer and 
product on cme side and fuel and product on the other. Their solution has 
instead of c^ + , and C^(“) instead of C^ + 

C^(-«) and C^(*) are the mass fractions of fuel and oxidizer in the volumes 
adjacent to the reaction zone. Following Gibson and Libby , ( 5 ) msy perhaps 
be interpreted as a mean rate of consumption of f lel at a point, related to 
the frequency of passage of reaction zones through that point. Equation ( 5 ) 
is also related to Spalding's models for premixed flames, in that w^ is 
proportional to e/k. The dependence on amount of reactants is different, 
however. It is evident that w^ as calculated by ( 5 ) does not go to zero 
when one of the reactants disappeacrs. For this reason w^ was set equal to 
zero when the sum of the mean value and the root -mean-square fluctuation of 
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mass fraction of either reactant vas less than iO « 

It 

The differential equations for and g^ also contain reaction rate 


terras, p and ^ * respectively. For the one-dimensional 

8 

shear flew studied by Chting , the species conservation, equation .for fuel can 
be vritten, after multiplying by and talcing the time average, as: 


Tlsi^ Chung's^ solution for the probability densities to evaluate vC^gives: 


so that ^ ^ equation for g^ by 

neglecting density fluctuations, giving p = p g^ . 


Summary of Eemations 

The foUc:’ing is a summary of the mean differential equations that 
were solved. They are vritten for axially symmetric flow in the boundary" 
layer approximation. Density fluctuations have been neglected. Various 
tutrbulence transport terms have been modeled as described by Launder and 
Spalding^*^. The conservation equations are: 


: 


&(p Ur)/Sx + 3(p VrySr = 0 


( 7 ) 


x-moimentum: 


p fj 3u/dx + p V 3u/3r = r“^(3/^r) (rjA^, 30/dr) - dp/dx (8) 



where 


)'i 


Energy: 


where 


Fuel: 


11 _ 

1 

8 ’ 

./• 

(9) 

p fj 3h^/3x + p f « r“^(3/3r)[|r(M-^ffjj)3E/5.r 

(10) 

+ xyiJS 


h=E + JS^+ic 

o 2 ' 


P U + P ^C^5r “ 

r~\W^r) [r(;i^/cp P . 

(u.) 


Beactioxx 
inyariant : 


p U P V 3C^/Sr “ r"\s/ar)[r(ii^/ff^)acyar] 


( 12 ) 


w'jtere 

. if. 


= ®0 - (VA'r*) 


Fuel 

fluctuations : 


P + P ^ = r“^(a/Sr)[r(iy'ag)3g^5rj (13) 


' °gl - Cg2 Sf/V^ * 2 p g^°*5 ^ 

The equation for oaygen fluctuations is the sane as (13) » hut with 
g^, C^, and replaced hy g^, C^> and ?P , respectively. Here, ^ 
■'equals (W^h/W^a) w^ . 


Turbulence kinetic 
energy: 


P U 3k/&x + p V 3k/3r = r"^(a/ar)|^r(^^Cj^) 3k/SrJ 

+ Hj(au/ar)“ - pt 


m 



9 


Turbulence _ - 

dissipation p 0 ftt/dx + p V 3e/3r ® r"^(3/3r) 3«/S2’J (l5) 

rate: 

P^(^U/ar)^ - C,2 P «^/k 

Tiiere is aI«o a differential equation for ^ of the same form as 
( 12 ). These differential equations vere supplemented by equations of state j 
which assuzae a mixture of perfect gases with variable specific heats. The 
numerical solution was performed using a finite Hiifference method based on 
that of Tataokeur and. Spalding . 

Ctomparisons with Experiment 

A limited evaluation of this reaction-rate model has been made by 

5 

comparing calculations with the experimental data of Kent and Bilger for 
velocity ratios U./U equal to 10, 2. These eo^riments w&re 

chosen as a standard of comparison because of the extensive flow field 
measurements, including initial conditions for u ♦ The experimental 
arrangement was a central hydrogen Jet Isstiing into a parallel, coaxial air 
stream. Both streams were subsonic and initially at about 3 OOK, There was 
a sli^t favorable pressure gradient, caused by the influence of the tuzmel 
walls . 

The initial conditions and the pressure gradient for the calculations 

5 

were the same as those measured . k and e were not measured directly, 
however. The initial values of k were (^termined by k — u^ /4 , 

l4 

based on conditions in the inertial sxiblayer of a boxmdaiy layer . Initial 
values of € were determined from k and the measxired U(y) through the 
approximation • “ k5u/9r. This was based on the appjroxamation tfmt 
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Reynold stress is proportional to k, and on (9) , initial conditions for 
vere assuBied to be g^ » 2k(C^0)^, and similarly for g^. 

The most extensive set of measurements reported by Kent and Bilger were 

for Ut/U = 10. In making calculations for this case, a number of values 

of the coefficient A in (5) were tried. The value 0.22 gave sotoewhat better ' 

comparison with experiment than the original estimate of O.17. Figures l->7 

show comparisons of calculations using A = 0.22 'with experiments for 

U^/U - 10 . Figure 1 shows the axial 'variation of centerline values of 

temperature and of mole fraction of Og, andiHgO. Figures 2-5 show 

radial profiles of these four quantities'' at x/D of ko, 80, 120, and I60. 

These comparisons show satl .factory agreement, overall, and are generally best 

nearer the injector (smaller x/o) an.i closer to the centerline (smaller v/tj) . 

The latter may be related to the omission from the analysis of Intermitt ency 

effects. Also, the omission of correlations in<rolv5.ng density fluctxiations 

is expected to be important at larger r/rj, where ther^ are large' density 

4 

gradients. Figures 3 and 5 also show resxilts calculated by Rhodes, et al , 
using an integral method, with equilibrium chemistry and their model for the ' 
mean density. At x/D = 80 (Figure 3), the present method and that of Rhodes, 
et al, are both In good agreement 'with experiment. At x/d = I60 (Figure 5)j 
however, the results of the present method are in substantially better 
agreement with e7qperime33t . 

Figure 6 shows flame contours. The and Og limits are where the 
mole fractions of Hg and Og are 0.01 and 0.005> respectively. The 
stoichiometric line is where and Og appear in stoxchiometrio proportions. 
If the intersection of the stoichiometric line with the axis is used as a 
measure of flams length, then the calculation oveiprcdicts this cose by about 
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'Ji», Turbulence levels on the axis ore shown in Figure 7. Bie calculated 
line is based on the assua^tion that equals (s/s) k, as for isotropic turbu- 
lence. 

Calculations were also made for Uj/0^ of 5 ^ind 2 using A equal to 

0.22. Figure 8 shows the axial yariation of temperature and of mole fractions 

of Og, and for Uj/U^ = 5- Eadial profiles of these quantities at 

x/d of 4o are shown in Figure 9* The seme information fbr Ut/u = 2 is shown 
' J 6 

in Figures 1C and 11, These comparisons show rather close agreement with 

experiment, as for - 1C. This lends some support to the idea that the 

coefficient A is nearly a constant, although this is not conclusive. Again 

there Is generally better agreement for smaller x/d and smaller r/rj. 

The axial variation of mean velocity for all thiee velocity ratios is 

compared with eaq^erlment in Figure 12. For U_/t)’ of 1C and calculations 

j e 

■ ' ■ 'i 

are within about of e:: 5 «rlment. Differences of up to l6^ occur for 

-2, however. This may be a result of the turbulence transport model, 

f.' 

because when the velocity excess in the jet is small, the model predicts a 

17 

velocity decay rate that is slower than measured . This applies for 

|j 

=2, where for exaiqple at x/d = 80, the calculation gave V“e = 1.052. 

// Although the reaction-rate model proposed here greatly simplifies the 

12,13 

actual physical behavior, this model and the turbulence transport model 

were generally successful in predicting measured values from the experiments 

5 

of Kent and Bilger . The favorable results obtained so far, plus the ease 
of use of the model, justify further consideration of this reaction-rate model 
for turbulent diffusion flfune calculations. 
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Suntnary 

The single reaction-rate moOel for turhuleirb diffusion flames px-oposed 
here has the following major characteristics; (l) the reaction rate is 
proportional to the tiu^bulence mixing rate} (2) the reaction rate is dependent 
on the mean mass fraction and on the mean square fluctuation of mass fraction 
of each reactant. 

Computations using this reaction-rate modrl have been generally 

successful in predicting meastired quantities for three hydrogen jet |b air 

5 

stream velocity ratios from the experiments of Kent and Bilger , 
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